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I. INTRODUCTION 



In recent years the study of the deeply virtual Compton scattering (DVCS) became one 
of the most popular topics in QCD due to the fact that it is determined by skewed parton 
distributions [1-3] which generalize usual parton densities introduced by Feynman. These 
new probes of the nucleon structure are accessible in exclusive processes such as DVCS and 
potentially they can give us more information than the traditional parton densities. In this 
paper we consider the small-x DVCS where the energy of the incoming virtual photon E 
is very large in comparison to its virtuality Q^. (The first study of the small-x DVCS was 
undertaken in Ref. [|]). To be specific, we calculate the DVCS amplitude in the region 

s > > -t > (1) 

where s = 2mE, m is the nucleon mass, and t is the momentum transfer. The DVCS in 
this region is a semihard processes which can be described by the BFKL (Balitsky-Fadin- 
Kuraev-Lipatov) pomeron 0. It turns out that at large momentum transfer the coupling 
of the BFKL pomeron to the nucleon is essentially equal to the Dirac form factor of the 
nucleon -Fi(t), so the DVCS amplitude in the region ([^) can be calculated without any 
model assumptions. The results obtained in this region can be used for the estimates of the 
amplitude at experimentally accessible energies where one or more conditions in Eq. (|Tp 
are relaxed. To be specific, we have in mind the HERA kinematics where x ~ 10^^ ^ 10"'^, 
> Q GeV^, and — t ~ 1 5 GeV^ ^j. Since there are only model predictions for the 
small-x DVCS in current literature , even the approximate calculations of the cross section 
in QCD are very timely. 

II. SMALL-X DVCS IN THE LOWEST ORDER IN PERTURBATION THEORY 

Similarly to the case of deep inelastic scattering (DIS), the amplitude of DVCS is deter- 
mined by the matrix element |p 
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= I dze^'^'^{p'\T{nz)fmp) (2) 

where q,p and q',p' are the initial and the final momenta of the photon and the nucleon, 
respectively. The momentum transfer is defined as r = p' — p. Since = —q^ is large 
we can use perturbation theory for the hard part of the DVCS process Q [l^. The typical 
diagram for the DVCS amplitude in the lowest order in perturbation theory is shown in Fig.l 
(recall that the diagrams with gluon exchanges dominate at high energies). It is convenient 




FIG. 1. A typical Feynman diagram for the high-energy 7*p — > scattering 
to calculate at first the imaginary part of the amplitude H"^^ 

V^"" = -ImT^^. (3) 

TT 

In the leading order in perturbation theory the amplitude at high energy is purely imaginary 
up to the ^ corrections (see e.g. the review At high orders in perturbation theory the 

amplitude will be purely imaginary in the leading logarithmic approximation (LLA) and we 
will restore the real part using the dispersion relations. 

At high energies it is convenient to use the Sudakov variables. Let us define the light-like 

2 

vectors pi = q', P2 = p' — ^Pi, then 

q = Pi{l - —) - xp2 - r_L q' = Pi 

s 
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p = P2{1 +x) -\ pi + r_L P = P2-\ Pi (4) 

s s 

where x = ^-j-^ — ^ — ^bj and t ~ —rj_ at large energies. Consider the integral over gluon 

momentum k = aupi + (3kP2 + 

where r + /c) and ($Ar)|^(fc, r + k) are the upper and the lower blocks of the diagram 

in Fig. ^ (stripped of the strong coupling constant g). Here a,b and ^,77 are the color 
and Lorentz indices, respectively. It is well known that in the Regge kinematics (= s ^ 




FIG. 2. Block structure of small-x DVCS in the leading order in perturbation theory 

everything else) ~ and ~ x so fc^ ~ —k'j_. Moreover, a's in the upper block are 
~ 1 so one can drop ak in the upper block. Similarly, /5's in the lower block are ~ 1 and 
one can neglect (3^ in the lower block. We get ($"^ = ^$^^): 

yAB^g^fd^^ ^ l^mk + r,-k) lm^fj'\-k-r,k) . (6) 
4nJ mn^kl{r + k)l ^"^^ ' ^ "^=0 ^ ^ ' ^ Pk=o 

At high energies, the metric tensor in the numerator of the Feynman-gauge gluon propagator 

reduces to g^'^ f so the integral @ for the imaginary part factorizes into a product 

of two "impact factors" integrated with two-dimensional propagators 
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where 



Mk^, ri.) = j-Apl I ^Im$«^%(-fc - r, k) 



ak=0 



(8) 
(9) 



and (j^efj is the sum of squared charges of active flavors {u,d,s, and possibly c). The 
photon impact factor is given by the two one-loop diagrams shown in Fig. R|. The standard 





FIG. 3. Photon impact factor 



calculation of these diagrams 12 1 yields 



where 



1 /"i da /"i da' ( i-, ^2 / 1 

/ — / \ Pfa'a' + Q^a'aa \ 

Jo 27r JO 27r J 



1 - 2«a)P|(e^, e^)± + 4aaa'[P^(e^, e^) - 2(e^, P)±(e^, P) J - 4aa(l - 2a)(r, e^)^(P, e 



(10) 



'111 



for the transverse polarizations A,B= 1,2 (cf. [|T^) and 
^ ' ^ 2g io 27r io 27r I ^ ^ J 



(12) 



(1 - 2aa)Pl{r, e^)± + Aaaa'[P'i{r, e")^ - 2(r, P)±(e^, P) J - 4aa(l - 2a)Q'\P, e")^ 



for the longitudinal polarization 



— {pi + XP2) 



(13) 
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. Here P± = k± + r±a and (a, b)± denotes the (positive) scalar product of transverse 
components of vectors a and b. At large transverse momenta ^ the impact factor 
(p!0|) reduces to 



/^^(A:^,r^)-^:4^^1n^. (14) 



47r2 Q2 ^ ^2 



The impact factor for the proton which decribes the pomeron-nucleon coupling cannot 
be calculated in the perturbation theory. However, in the next section we demonstrate that 
at high momenta k'j_ ^ this impact factor reduces to 

/iv(fc±,r^)''="^Vr"(t) (15) 

where Ff'^"(t) is the sum of the proton and neutron Dirac form factors. As we shall see 
below, the characteristic transverse momenta in our gluon loop are large so the estimate 
(P^) is sufficient for our purposes. Substituting the nucleon impact factor (|T^) into Eq. ([Zp 
we obtain 

V l^A:i(r + fc)i- ^^^^ 

Performing the final integration over k^, one gets 

V^"" = 

flavors 



X \ TV 

((e^, e^)x In' ^ + ^ j " """^^ + i^'^' ""^^ + ^^^/^'^ ) ^^^^ 

for the transverse polarizations and 



3B 



^ ^ ^ ^ flavors V 1*^1 |r| / d J 



for the longitudinal one. The longitudinal amplitude (p^q) is twist-suppressed as ^ in 



Q 

comparison to the transverse amplitude (|T^ (as it should, due to the fact that t ^ 
corresponds to real incoming photon). 
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Since the integral over k± ([T6|) converges at fc_L ~ Q the region k± ~ m, where we do not 
know the nucleon impact factor, contributes to the terms ~ 0(t/Q^) which we neglect. 

III. NUCLEON IMPACT FACTOR 

In the lowest order in perturbation theory there is no difference between the diagrams 
for the nucleon impact factor shown in Fig. |^ and similar diagrams with two gluons replaced 



r+k k 




FIG. 4. Nucleon impact factor 

by two photons (up to the trivial numerical factor cp = ^ and replacement oi e ^ g). In 
this case the lower part of the diagram can be formally written as follows: 

^^{-k - r, k) = -r,k) = '^-tp^.p'i J dze^'^p'\T*{J,{z)MO)}\p) (19) 

where = u'~ffj,u + d'-f^^d. The T* means the T-product where the diagrams with pure gluon 
exchanges in t-channel are excluded; by definition, such diagrams contribute to subsequent 
ranks of BFKL ladder rather than to impact factor. (This is the reason why we have not 
included in J the contribution of strange quarks). Since A;^ in our case is large and negative 
{-k"^ = k\'^ im?) we can expand the T-product of two currents near the light cone (see e.g. 

0) 

^M{k,r + k) = |- / dze'^^^^{p'\-i,{z)[z,Q] ^iV^(O) + V^(0)[0, ^] M{z)\p)l,^, (20) 

where again (...)* stands for the matrix element with pure gluon exchanges ex- 
cluded. Here [x^y] denotes the gauge hnk connecting the points x and y ([x, y] = 



7 



Pexp (ig Jq du{x — yYA^{ux + (1 — u)y^). The matrix element (16) can be parametrized 
in terms of skewed parton distributions as follows 



{p\\'\q{z)[z,Q] ^ig(0)b,A):.=o= (21) 

(p', A') M{P, A) C dXe'^^-'^^P'V^.iX, t) + -^uip', A') f^p, A) dXe'^^-^'^^'WHX, t) 
Jo 2m Jo 

{p',X'\m[0,z] ^iq{z)\p,)^):^=o = 

ip', A') MiP, A) /' dXe-'^'P^V'^iX, t) + ^u{p', A') /Mp. A) /' dXe-'^'P'W^^iX, t), 
JO 2m Jo 



u 



u 



where V^(X, t) and W"(X, t) are the nonflip and spin- flip skewed parton distributions for 
the valence u quark (recall that we must take into account only valence quarks since we 
forbid diagrams with pure gluon exchanges). Similarly, V^(X, t) and W^(X, t) refer to the 
valence d-quark distributions. At large energies u{p' , A') ;^iu{p, A) = s6x\i, so 



{p',X'\m[0,z] ^M^)-q{z)[z,0] ;^ig(0)b,A):.=o 



(22) 



g -iXpz _ ^i(X-x)pz 



1 



/.«( 

After integration over z the lower block (p!9| ) reduces to 



s5aa'V,^(X, t) + —u{p\ A') ^1 fMP. A) t) 
2m 



^N{-k - r, k) 



(23) 



2 /-i 
-/ rfX 

3s JO 



(X - x)s + 2pi-k -Xs + 2pi-k 

-P - 2p ■ k{X -x)-ie ~ -k^ + 2p ■ kX - ie 



6,y{v:{x,t) + v',{x,t)) + ^Hp'A') Af±u{p,X){W:{X,t) + W',{X,t)) ] . 



The nucleon impact factor (^ is the integral of the imaginary part of r.h.s. of eq. (^3]) over 
energy 

lN{k±, r_L) = / — — Im$Ar(-(afc )pi - k^ - r±, atPi + fc±) = (24) 

Jo ^TT S 

- / dak I dX s{X — x)5{k']_ — a^s^X — x)) — sXdik^ + aksX) 
3 Jo Jx 

6xx'iV:iX, t) + Vi{X, t)) + -^u{p\ A') A f±uip, A)(W:(X, t) + >V,^(X, t)) 

2ms 

- ^ £ dX (5aa'(V:(X, t) + Vt{X, t)) + ^^(P, A') A Mp, A)(W:(X, t) + W'^iX, t))^ 
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Since valence quark distributions decrease at a; — we can extend the lower limit of inte- 
gration in r.h.s. of eq. (|2^ ) to and obtain 

/7v(fc±,rx) (25) 
i J^' dX {dMmX, t) + (X, t)) + i^u{p\ A') ^1 fMP, A)(>V:(X, t) + Wt{X, t))) . 

Let us recall the sum rules ^ 

f dX {Tl{X, t) - Tl{X, t)) = F?{t) 
Jo 

f dX {ICliX, t) - /Cf (X, t)) = Flit) (26) 

where J^^(X, t) and /C^(X, t) are the total (valence + sea) nonfiip and spin- flip skewed quark 
distributions while J^|(X, t) and /Cf (X, t) are the antiquark ones. Here Flit) and F^it) stand 
for the g-quark components of the Dirac and Pauli form factors of the proton). Since the 
contribution of sea quarks drops from the difference JF'' — JF^ we can rewrite eqs. p6| ) as the 
sum rules for valence quark distributions 

/' dXVl{X, t) = F?{t), f dXWl{X, t) = Flit). (27) 

Jo Jo 



Substituting this estimate to eq. ( pSf ) and using the isospin invariance, we get the final result 
for the nucleon impact factor at large transverse momenta 

iNik^, rx) Sxx'Fr'^it) + -^uip', A') /Mp, A)F|+"(t), (28) 

2ms 

where Ff+"(t) = Ff (t) + F{^(t) and F|+"(t) = F|(t) + F^it). As usual, Ff^^^ and F^^"^ are 
the Dirac and Pauli form factors of the proton (neutron), respectively. With our accuracy 
they can be approximated by the dipole formulas 

+ ^Fi= G% = . Ff + Ff = G% = -]l . 

I "^0.7GoV2 ) l^0.71GcV2 J ^29) 

Fn I t rrn /^n n T?n- i I?n- —1.91 
1 + - - U i^l + - ^AI - 7- ^1 ^ 

l^+0.71GevSJ 

which leads to Q 



*Literally, one obtains 
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= rV^' ^'''^ = ° ^^^^ 



^0.7GeV^, 

. Note that the spin-flip term turned out to be neghgible for our values of t. Moreover, it 
vanishes at t = which suggests that it is numerically small at all t. 
Our final estimate of the nucleon impact factor is 

lN{k^,r^)''^="'' hyFnt) (32) 

where Ff"*"" is given by the dipole formula (^) f\. In what follows we shall omit the factor 
6 XX' (as it was done in eq. (|15|)) since all our amplitudes will always be diagonal in the 
proton's spin. 



j,p+n. . _ 1 1 + 0-88 jg^ _ 0-12 . . 

( \t\ \2 1+ 1*1 ' ^2 - , x2 l^^i 



but with our accuracy we can use the estimate (pll). 

^ The dipole formula for the neutron form factor does not seem to work as well as the dipole 
formula for the proton form factor. As a measure of the uncertainty we can compare the results 
obtained from eq. ( plj ) to those obtained using the model from Ref. (which was fit only to the 
proton form factor) 



<"(i) = \ I' dX {V^{X, t) + V'AX, t] 
V^{X, t) = imX~^-^X^-^{l + 6X) exp ' ^ 



X2.8GeV2y 

Vi{X,t) = 0.54X-0-6x4-2(l + 8X)exp (^-|_^^ (33) 

The results for the DVCS cross section in this model are about 1.5 times bigger than the results 
obtained from the dipole formula (31). 
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IV. THE BFKL LADDER 



In the next order in perturbation theory the most important diagrams are those of the 
type shown in Fig. . Calculation of this diagrams in the leading log approximation 



k+r 



k'+r 




FIG. 5. Typical diagram in the next-to-leading order in perturbation theory 



yields 



AB 



V 



^(Ee?)(6«.lnl)/. 



47r2 47r2 P^{r + k)\ 
where K{k±, k'j_,r±) is the BFKL kernel 



(A;^)2(r + k')l 



(34) 



K{k^,k'^,r^] 



2 1, 



(fc-fc')i 



(k-k') 



-k'±)Idp^{^;r^ + j;;^^^i^] (35) 



(p-r)Uk-py. 



As we shall see below, the integral over k'j_ converges at \k'j_\ ^ m so we can again use the 
approximation (^) for the nucleon impact factor. One obtains 



•f- Actually, this diagram gives the total contribution in LLA if one replaces the three-gluon vertex 



in Fig. (|5|) by the effective Lipatov's vertex |11] 
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and therefore the amphtude (|34D takes the form 

TT 



{k - r)^ 



(36) 



Finally, the integration over k yields 

^ flavors ^ ^ ^- 

^^"^^ ^^ln^|| + 21n||-2 + C(3)] + (^(e^r)^(e«,r)^-(e^e-)^) ) (38) 



2 /a 



where the accuracy is 0{j^). 

In the next order in BFKL approximation (see Fig. ^) it is still possible to obtain the 
DVCS amplitude (^ in the explicit form (we have not obtained the explicit expressions for 
higher-order terms in the BFKL expansion ( ^8|) |). The amplitude in this order is 

V'--^(i:e?)(6".l.>i)7=g^^^/(*x.rJ ^ (39) 

Once again, if we use the fact that the integral over k'j_ converges at ^ m we can 

approximate the nucleon impact factor by eq. (^), and obtain 
c/2K r d'^k" I 1 



■-K{k^, r^) ^ i^(fc"^, fcl, r^) ^ /^(fcl, rl) 



47r2 7 47r2 ^ (fc")^(r + fc")i ' ^' (A;')i(r + A;' 

1 pp+n... f d'k^\Kik^,k^\,r^) ( ik'\y (fc"- 
47r^ ^ 47r2 (A;")^(r + A;")i V ^ ri 



Ff+"(t) In^ ^ + In^ ' . (40) 



167r2 \ 



^It is possible to write down the result of the summation of the BFKL ladder in the form of Mellin 
integral over complex momenta using the Lipatov's conformal eigenfunctions of the BFKL equation 
in the coordinate space. Unfortunately, we were not able to perform explicitly the integration of 
the Lipatov's eigenfunctions with impact factors and without it the Mellin representation of the 
DVCS amplitude is useless for practical applications. 
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FIG. 6. Typical diagram in the next-to-next-to-leading order in perturbation theory 
The resulting integration over k± yie Ids 




As we mentioned, we were not able to obtain the explicit expressions for the amplitude in 
higher orders in perturbation theory. It turns out, however, that for HERA energies the 
achieved accuracy is reasonably good; the estimation of the next term gives ~ 30% of the 
answer at not too low x (see the discussion in next section). Our final result for the DVCS 
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amplitude with transversely polarized photons is Q 

= (42) 



2 fa,iQ)Y 



X \ 71 

where 



flavors 



\ r1 



2 

TT X ' 2 \ TT X 



.'(x,QVt) = l + ^^ln- + U^^ln-l • (44) 



Note that the spin- dependent part ~ v' does not contain any In ^ and is hence much smaller 
than the spin-independent part ~ v. For the longitudinal polarization ([T3| ) the amplitude 
is twist-suppressed as ~ so we have not calculated any terms beyond eq. (p!8|). In the 
numerical analysis carried out in next sections we keep only the spin-independent part of 
the amplitude 

^ \T.eielV^' = I (^)'( E el)Fr{t)vM\t). (45) 

\ / flavors 

The above expressions give us the imaginary part of the DVCS amplitude. For the 
calculation of the DVCS cross section we need to know also the real part of this amplitude 
which can be estimated via the dispersion relation. For the positive-signature amplitude if^ 
{=\Y.ei_elH^^) we get p[ (see also §) 

ReHAs) = -tiin\s^\\mHAs), (46) 
2 \ as J 

which amounts to the substitution 



**In the leading logarithmic approximation it is not possible to distinguish between as{Q) and 
as(-y|t[) - to this end one needs to use the NLO BFKL approximation [l^ (see also which is 
beyond the scope of this paper. 
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Ins ^ -( ln(— s — ie) + Ins) 



(47) 



in our amplitude (|i5|). Thus, the real part is 



R = -ReHj 

TT 

r{x,Q',t) = 



2 

X \ TT 



(E el){Fm + Fmr{x,Q\t) 



flavors 



TT 



rr \Q t t ' 



(48) 



V. COMPARISON WITH THE DEEP INELASTIC SCATTERING 

It is instructive to compare the DVCS amplitude V^^ given by eq. with the corre- 
sponding amplitude for the forward 7* scattering 

T^^ = teX J dze''^^{p\T{f{z)f{0)}\p). (49) 

The imaginary part of this amplitude is the total cross section for deep inelastic scattering 
(DIS) 



TT 



pApB 



(50) 



For example W averaged over the transverse polarizations of the photons is 

W± "^l^^^fW^'' = F^i^^Q') = ^^2(0:, Q^) (51) 

(at the leading twist level we have the Callan-Gross relation F2 = 2xFi). We will compare 
the imaginary part of the DVCS amplitude V± given by eq. (^51) to the result for W± 
calculated with the same accuracy. (We use the notation Wj_{x) rather than Fi{x) in order 
to avoid confusion with Fi{t)). 

Similarly to the DVCS case, the DIS amplitude has the form (cf. eqs.(|16|,)(p3), and (|39D): 
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(pk^ 1 



/±(A;±,0) 



1 + 14 In - / d^k'^K{k^, k'^, 0) + 

In'^Jd^k'^ I d''k'\K{k^,k'\,0)-^K{k'\,k'^,0) 



1287r 



r^lNik'^,0) (52) 



where I±{k±, 0) is the virtual photon impact factor averaged over the transverse polarizations 



Hi 



1 da /-i rfa' kl{l - 2aa)(l - 2aV) 



(53) 



2 7o 27r 7o 2tt kj_a'a' + Q'^a'aa 

The nucleon impact factor lN{k'j_,0) cannot be calculated in perturbation theory since it 
is determined by the large-scale nucleon dynamics. However, we know the asymptotics at 
large k± ^ m 



lN{k±,0) 



Ff+"(0) = 1 



(54) 



Also, at lN{k±, 0) ^ at A; ^ due to the gauge invariance. It seems reasonable to model 
this impact factor by the simple formula 



lN{k±,0) 



(55) 



which has the correct behavior both at large and small k±. With this model, the DIS 
amplitude ( ^21) takes the form 



F2 

W^ = — = 
2x 

4 fasiQY 

3X I TT 



(56) 



E 



flavors 

3«.i„l ri,.3Q^ , 7 ,^,Q^ , 77,^0^ , 131 

131 



2 m"* 6 m"* 18 / 



In- -ln^:^ + -ln^:^ + -ln:^ + — + 2C(3) 
TT a; \ o m"^ 12 m"^ 18 27 



1 /3a,, , 1 
m — 



2 \ TT X 



24"' m2 ' 36"^ m2 ' 36"^ m2 + ^ 27 + ^^^^^^ ^2 + 



g2 1396 vr^ 14 



15 + T<P' 



Note that the coefficients in front of leading logs of Q^, determined by the anomalous 
dimensions of twist-2 operators, coincide up to the factor 2/3. The graph of the model (|56D 
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F2 



Q**2=10 



F2 



Q**2=35 




Log X 




-3 -2 



Log X 



FIG. 7. F2{x) from eq. ( |56| ) versus experimental data at Q'^ = lOGeV^ and Q"^ = 35GeV^ 

versus the experimental data is presented in Fig. 7 for = lOGeV^ and = 35GeV^ (we 
take Ee^ = f ). 

In the case of DIS it is possible to calculate explicitly the next term in BFKL series ( ^Gf) 
f]\ . It has the form 



4 fasiQY 

3X \ IT 

+3C(3))ln^ 



(E ' 

flavors 



+ 



+ 



6 V TT xj \120 m2 144 10 



,131 



54 



+ 



,1396 TT^ 



+ 7C(3))ln:^ + 



4736 77r^ 



77 



81 15 " m? 243 90 3 

\3 



C(3) + 6C(5) 



(57) 



The ratio of this (a^ In x) term to the sum of the first three ones ( p6|) is presented in Fig. ^ 
for = lOGeV^ and = 35GeV^. From these graphs we see that the sum of the first tree 
terms gives the reliable estimate of the DIS amplitude at not too low x and it is expected 
that the same will also be true for DVCS amplitude 0. 



^For DIS it is possible to write down the total BFKL sum as a Mellin integral and unlike DVCS 
the integrals of impact factors with the BFKL eigenfunctions (/c^,)" 2+*'^ can be calculated explicitly. 
Eqs. (^) and (57) correspond to the expansion of this explicit expression in powers of aslnx. 

At very small x ~ 10^^ 10~^ the full BFKL result for F2 in our model is growing more rapidly 
than Fig. |^. On the other hand if one takes into account the NLO BFKL corrections [16| the 
result for F2 at very small x goes well under the experimental points. This indicates that at such 
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0.4 



Q**2=10 



-2.5 -2 -1.5 -1 



Log X 



Q"2=35 




-2.5 



-1.5 



Log X 



FIG. 8. The ratio of the In^ x term to eq. (|6D at = lOGeV^ and = ssQeV^ 



It is instructive to compare the t-dependence of our DVCS amplitude P3| ) with the model 
used in the paper [0] 



1 



^2(x,t,g^) = iFi(x,g^) ^ 

V 0.71^ 



(58) 
(59) 



where R ~ 0.5 for our energies. (Literally, the model used in ref. corresponds to Vi 
but it is more natural to approximate the t - dependence by the dipole formula ||2^). The 
comparison is shown in Fig. for = lOGeV^, Q"^ = 35GeV^ and x=0.01, x=0.001. 



VI. DVCS CROSS SECTION 



In order to estimate the cross section for DVCS at HERA kinematics (Q^ > 6GeV^ and 
X < 10"^) we will use formulas from Ref. 0] (see also Ref. |^) with the trivial substitution 
^F2(x)i?~^e^*/^ — > V±{x,Q'^,t). The expressions for the DVCS cross section, the QED 



X we need to unitarize the BFKL pomeron, which is currently an unsolved problem. (The best 
hope is to find the effective action for the BFKL pomeron (see e.g. |2C|, [^)). On the contrary, at 
"intermediate" x ~ 0.1 -r- 0.001, we see from Fig. |^ that, since the corrections almost cancel each 
other, it makes sense to take into account only a few first terms in BFKL series. 




x=0.01,Q**2=35 



x=0.001,Q**2=35 





FIG. 9. The ratio Vi/V^ (lower curve) and V2/Vj_ (upper curve). 
Compton (Bethe-Heitler) cross section, and the interference term have the form [y = 1 —y) 



dxdydtd(j)r 

^^QEDC 



vra X- 



1+y 



y{l + y"^ 



dxdydtd(j)r tcx \t\Q'^y 



{Vl{x,Q\t) + Rl{x,Q^t)) 



dxdydtd(j)r 



Q\y\t\ 



Am? 



R±{x,Q^,t)Ff{t) cos0^,. 



(60) 
(61) 
(62) 



Here y = 1 - ^ {E and E' are the incident and scattered electron energies, respectively, 
as defined in the proton rest frame) and (pr = <pe + <Pn where is the azimuthal angle 
between the plane defined by 7* and the final state proton and the x — z plane and 0e 
is the azimuthal angle between the plane defined by the initial and final state electron and 



"The expression for the interference term from ref. is corrected by factor 2 |22], |25| 
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FIG. 10. Asymmetry versus y = 0.1 0.6 and |t| = 1-^5 GeV^. 

At first let us discuss the relative weight of the above cross sections. We start with the 
asymmetry defined in ref. 124 



C rf0.rfaDQi ^^^^ 

where 

da^Qi = rfa^vcs + ^^QEdc ^ ^^int^ (g4) 

The asymmetry shows the relative importance of the interference term, which is proportional 
to the real part of the DVCS amplitude. In our approximation the asymmetry is 
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4vr2(Ee2)2(t;2 + r2)(^) (l + 



4m? J y\t\ 



(l + 7.84 



(65) 



The plots of asymmetry versus y and \t\ are given by Fig. |T0 



Second, we define the ratio of the DVCS and Bethe-Heitler cross sections [0 



Diy,t) 



dac 



(66) 



'QEDC 2/2 (l + 7.84^) 

This ratio is presented on Fig. 11. We see that there is a sharp dependence on y; at y > 0.2 



x=0. 01,0**2=10 




x=0.001,Q**2=10 




FIG. 11. The ratio D{x, Q'^/t) versus y = 0.1 ^ 0.6 and |t| = 1 ^ 5 GeV^ 

the DVCS part is neghgible in comparison to Bethe-Heitler background whereas at y < 0.05 
the QEDC background is small in comparison to DVCS. 
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Finally let us estimate the relative weight of the DVCS signal (starting from |t| = 1 
GeV^) as compared to the DIS background. We define (cf. ref. ^) 



A-Ka 



0"(7* + p ^ 7* + p) 



(67) 



At g2 = lOGeV^ we find = 1.56 x 10"^ for x = 0.01 and R^ = 2.36 x 10"^ for x = 0.001, 
while for = SSGeV^ we find R^ = 0.62 x 10"^ for x = 0.01 and R^ = 0.71 x 10"^ for 



X = 0.001. 



Asymmetry at x=0 . 01, Q**2=10 




Asymmetry at x=0 . 001 , Q* *2=10 




FIG. 12. Asymmetry with the correction factor (68). 

The expressions (60)-(62) are correct if <^ \t\ up to 0{^) accuracy with the notable 
exception of the correction 0(-^^) coming from the expansion of electron propagator in the 
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u-channel of the Bethe-Heitler amplitude. As suggested in ref. at intermediate t one 
can keep the propagator in unexpanded form (and expand the rest of the amphtude, as we 
have done above). This amounts to the replacement 



y^y 



HI ^ , I^IVn (2 — y) / \t\ , ,\t\ r, 



(68) 



in the numerator in eqs. (61) and (62) (see ref. [^). The resulting asymmetry (63) is 



x=0 . 01, Q**2=10 




x=0. 001,0**2=10 




FIG. 13. The ratio D{x,Q'^ /t) with the correction factor (68). 

presented in Fig. 12. We see that the correction factor (68) crucially changes the behavior 
of the asymmetry due to the fact that it restores the azimuthal dependence of the QEDC 
amplitude which was not taken into account in eqs. (60-62). In order to find asymmerty at 
these and t with greater accuracy one should take into account other twist-4 contributions 
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as well. On the contrary, the ratio D{x,Q'^ /t) does not change much (see Fig. 13) so we 
hope that our leading-twist results for the ratio presented in Fig. 11 are reliable. 



VII. CONCLUSION 

The DVCS in the kinematical region is probably the best place to test the momentum 
transfer dependence of the BFKL pomeron. Without this dependence, the model (^) would 
be exact, hence the upper curves in Fig. ^ indicate how important is the t-dynamics of the 
pomeron. We see that the t-dependence of the BFKL pomeron changes the cross section at 
t > 2GeV^ by orders of magnitude and therefore it should be be possible to detect it. 

The pQCD calculation of the DVCS amplitude in the region (|I|) is in a sense even more 
reliable than the calculation of usual DIS amplitudes since it does not rely on the models for 
nucleon parton distributions . Indeed, all the non-perturbative nucleon input is contained 



in the Dirac form factor of the nucleon [***] , which is known to a pretty good accuracy. (Of 
course any reasonable models of nucleon parton distributions such as ( p^ should reproduce 
the form factor after integration over X). 

Finally, let us discuss uncertainties in our approximation and possible ways to improve it. 
One obvious improvement would be to calculate (at least numerically) the next ~ (a^ Inx)^ 
term in the BFKL series for the DVCS amplitude. Hopefully, it will be as small as the cor- 
responding calculation of the DIS amplitude suggests. Second, there are non-perturbative 
corrections to the BFKL pomeron which we mention above. These non-perturbative correc- 
tions correspond to the situation like the "aligned jet model" when one of the two gluons in 
Fig. (P is soft and all the momentum transfers through the other gluon. It is not clear how 
to take these corrections into account, but one should expect them to be smaller than the 
corresponding corrections to F2{x) which come from two non-perturbative gluons in Fig. 1 



*** There are, of course, the non-perturbative corrections to the BFKL pomeron itself. At present, 
it is not clear how to take them into account. 
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(in other words, from the "soft pomeron" contribution to ^2(2;)). 

The biggest uncertainty in our calculation is the argument of coupling constant which 
we take to be . As we mention above, it is not possible to fix the argument of as in the 
LLA, so we could have used as(|t|) instead. We hope to overcome this difficulty by using 
the results of NLO BFKL in our future work. 
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